Speed of disentanglement in multiqubit systems under depolarizing channel * 
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We investigate the speed of disentanglement in the multiqubit systems under the local depolar- 
izing channel, in which each qubit is independently coupled to the environment. We focus on the 
bipartition entanglement between one qubit and the remaining qubits constituting the system, which 
is measured by the negativity. For the two-qubit system, the speed for the pure state completely de- 
pends on its entanglement. The upper and lower bounds of the speed for arbitrary two-qubit states, 
and the necessary conditions for a state achieving them, are obtained. For the three-qubit system, 
we study the speed for pure states, whose entanglement properties can be completely described by 
five local-unitary-transformation invariants. An analytical expression of the relation between the 
speed and the invariants is derived. The speed is enhanced by the the three-tangle which is the 
entanglement among the three qubits, but reduced by the the two-qubit correlations outside of the 
concurrence. The decay of the negativity can be restrained by the other two negativity with the 
coequal sense. The unbalance between two qubits can reduce speed of disentanglement of the re- 
maining qubit in the system, even can retrieve the entanglement partially. For the fc-qubit systems 
in an arbitrary superposition of Greenberger-Horne-Zeilinger state and W state, the speed depends 
almost entirely on the amount of the negativity when k increases to five or six. An alternative 
quantitative definition for the robustness of entanglement is presented based on the speed of dis- 
entanglement, with comparison to the widely studied robustness measured by the critical amount 
of noise parameter where the entanglement vanishes. In the limit of large number of particles, the 
alternative robustness of the Greenberger-Horne-Zeilinger-type states is inversely proportional to k, 
and the one of the W states approaches 1/yk. 

PACS numbers: 03.67.Mn, 03.65.Ud, 03.65.Yz 



I. INTRODUCTION 

Entanglement [1], which reveals essentially difference 
between quantum and classical mechanics, plays a vital 
role in many discussions on the fundamental issues of 
quantum mechanics. In recent years, it has been found as 
a resource of quantum communication and computation 
[2], such as quantum key distribution [3, 4], quantum 
teleportation [5] and quantum dense coding [6]. 

However, it's impossible that a quantum system iso- 
lates from its environment completely, so the system un- 
avoidably loses its coherence due to interactions with the 
environment. In composite systems, entanglement is an 
manifestation of quantum coherence, which will decay 
under the influence of decoherence. This dynamics of en- 
tanglement in the open systems has recently attracted 
the interest of many researchers [7-19]. It has been 
shown that in some cases entanglement can vanish in 
finite times. This phenomenon of sudden loss of entan- 
glement is known as entanglement sudden death (ESD) 
[7-11, 20]. The experimental evidences of the interesting 
phenomenon have been reported for atomic ensembles 
[12] and optical setups [13, 14]. 

One of the topics in this direction is about the ro- 
bustness of entanglement, in which the stability of an 
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entangled system is a major concern to the researchers. 
Recently, many investigations on the robustness, espe- 
cially of the multiqubit quantum systems, have been 
reported. The concept for the robustness of entangle- 
ment has been introduced by Vidal and Tarrach [21] as 
a measure of entanglement corresponds to the minimal 
amount of mixing with separable states which washed 
out all entanglement. Then, in some articles [15-19], the 
robustness is defined as the critical amount of decoher- 
ence, where the entanglement entanglement vanishes. In 
[15], the authors pointed out that, under depolarizing 
channels, the robustness of entanglement for multiqubit 
systems increased with the qubit number k. Aolita et 
al. [16] showed that the time, when the entanglement 
of a pure Greenberger-Horne-Zeilinger (GHZ) state sys- 
tem became arbitrarily small and useless as a resource, 
was inversely proportional to the particle number fc, al- 
though its ESD time increased with k. In [17], Man et 
al. studied the entanglement robustness measured by the 
time when an entangled state became completely sepa- 
rable or remained entangled with a negligible entangled 
amount, and showed that entanglement robustness can 
be enhanced by local operation for the generalized k- 
qubit GHZ-type states. Zhao and Deng [19] studied the 
multiqubit systems under local depolarizing channel, and 
found a residual effect on the robustness of a three-qubit 
system in a superposition of GHZ state and W state. 

Inspired by the results given by Zhao and Deng [19], 
the main purpose of this paper is to study the relation be- 
tween the stability of an entangled system and its entan- 
glement properties. However, the robustness, measured 
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by the critical noise parameters where the ESD occurs, 
is determined by the nature of the series of intermediate 
states from being entangled to separable under the pro- 
cess of decoherence. In addition, as shown in [16], the 
ESD time is irrelevant to assess the robustness of mul- 
tiparticle entanglement in some channels. Therefore, we 
present the linear response of the entanglement under an 
infinitesimal noise, namely the speed of disentanglement 
(SoDE) , as a quantitative signature of the resistance of an 
entangled state. On the other hand, the behavior of the 
entanglement at the beginning of the decoherence process 
becomes very importance in the application for quantum 
information, since the coherence time of a qubit has been 
greatly prolonged in some recent schemes [22, 23]. It is 
not the first time to concern about the time-derivative 
of the entanglement, for instance, the derivative of the 
tangle is studied in the investigation about the entangle- 
ment flow in multipartite systems [24]. To exclude the 
properties of an entangled system affecting the SoDE be- 
side the entanglement, we adopt the model of a multi- 
qubit system in which each qubit is coupled to its own 
depolarizing environment individually, the same as that 
in [15, 19]. The model is invariable under both the local 
unitary (LU) transformations and the permutations of 
the qubits. Furthermore, in the limit of large number of 
particles, the qualitative behavior of an entangled state 
is largely independent of the specific decoherence model 
[25]. 

As the first trial, in the present paper, we only consider 
speed of the bipartition entanglement between one qubit 
and the remaining k — 1 qubits constituting the fc-qubit 
system. We adopt the negativity N [26-28] to measure 
the entanglement, which is widely used in recent investi- 
gations of the entanglement dynamics in multiqubit sys- 
tem [15, 19, 29, 30]. Compared to other entanglement 
measurements, such as the concurrence C [31, 32], a dis- 
tinct advantage of the negativity is its computability. Be- 
sides, for the arbitrary bipartite pure states, N > is the 
necessary and sufficient condition for entanglement [32]. 
Most of the recent works about the robustness of entan- 
glement in multiparticle system focus on the pure states 
case. It is certainly true that our definition of the SoDE 
is not restricted to the negativity. However, we find the 
SoDE for concurrence are very close to the one for neg- 
ativity, in some special cases where the concurrence or 
its bound can be calculated. The details are given in the 
conclusion and discussion section. 

This paper is organized as follows. In Sec. II we make a 
brief review for the the decoherence model and the entan- 
glement measurements, and give some general formulas 
for the SoDE. In Sec. Ill and IV, we study the SoDE for 
two-qubit arbitrary states and three-qubit pure states, 
respectively. In Sec. V, we give the result for several 
special states in multiqubit system, such as GHZ-type 
state, W-type state, etc. Conclusion and discussion are 
made in the last section. 



II. DEFINITIONS AND GENERAL FORMULAS 

A. Decoherence Model 

Under local noise environments, without the interac- 
tion between the subsystems in the fc-party system, the 
dynamics of each particle is governed by a master equa- 
tion which depends on its own environment [29]. From 
the master equation, one can obtain a completely positive 
trace-preserving map £j which describes the evolution of 
the corresponding subsystem [29]: pi(t) — £^(0), (i = 
l,...,fc), and for the whole state p(t) = <8'J : =1 eip(0). In 
the Born-Markovian approximation, the channel can be 
described by its Kraus representation [19, 29, 33] as 

M 

where Eji, j = 0, M are the so-called Kraus operators 
needed to completely characterize the channel. 

In detail, for the one-qubit quantum system, the Kraus 
operators of the partially depolarizing channel can be 
expressed as 



E = y/1 -p>I, Ej = 




where I is the unit matrix, p' = with p being the 
depolarization parameter, and Oj (j = 1,2,3) are the 
corresponding Pauli matrices. In the Bloch sphere repre- 
sentation for the qubit density operator, the map can be 
written as ep — \{I -\- sr ■ <r), where s = 1 — p. Following 
the Ref. [34], we consider the depolarization parameter 
s = s(t) = e~ Kt , where k is a decay constant determined 
by the strength of the coupling to the environment and 
t is the interaction time. Without loss of generality, we 
set k = 1 and s — e~* in the present work. 



B. Entanglement measurements and invariants 

Several concepts have been presented to quantify or de- 
scribe the entanglement in quantum systems, such as the 
entanglement of formation [31], the entanglement cost 
[35], the distillable entanglement [36], the relative en- 
tropy of entanglement [37], and so on. Here, we briefly 
review the concurrence [31] and negativity [26] for bipar- 
tition entanglement, and the entanglement invariants for 
three-qubit pure states, which are utilized in this paper. 

The negativity as a measure of entanglement was intro- 
duced by Zyczkowski et al. [28] , base on positive partial 
transpose (PPT) criterion [2 7], which is necessary for sep- 
arability, but sufficient only for 2 ® 2 and 2 ® 3 systems 
[38] . For a bipartite system described by the density ma- 
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trix p, the negativity is defined as [27, 28] 



(3) 



where Xj are the negative eigenvalues of p T and T denotes 
the partial transpose operation on one of the subsystems. 
For a pure m <8> n (m < n) state in the standard Schmidt 
form 



= y^ j a%\lJ.iVi) 



(4) 



where a» S [0, 1] (i = 1, m) are the Schmidt coefficients 
satisfying ^ of = 1, [ft) and are the orthonormal 
basis of the two subsystems, the negativity is given by 
[32] 



mm 



(5) 



from which one can find the fact that, Af > is the 
necessary and sufficient condition for the entanglement 
in pure states. 

The concurrence originated in the investigation of the 
entanglement of formation [31]. The concurrence of a 
pure state \tp) in a bipartite system is given by [32] 



i 



Trpi) 



\/2( 



1 



Trp|), 



(6) 



where qa = TtbIV'XV'I is the partial trace of \ip)(ip\ over 
subsystem B, and has a similar meaning. For a mixed 
state, the concurrence is defined as the average concur- 
rence of the pure states of the decomposition, minimized 
over all decompositions of p = ^2jPj\i^j){i>j\, 



C(p)=mmJ2PjC(\^)). 



(7) 



For the two-qubit case, it is equivalent to the entangle- 
ment of formation [31] and can be expressed explicitly 
as 



C(p) = max{0, Ai - A 2 - A 3 - A 4 }, 



(8) 



in which Ai, A4 are the square roots of the eigenvalues 
of the operator R = p(a y ® o y )p*{o y <g> a y ) in decreasing 
order and o~ y is the second Pauli matrix. 

For an m (£> n (m < n) arbitrary quantum state p, 
the negativity is proved to provide a lower bound of the 
concurrence as [32] 



C( P ) > J— -Af(p). 

y m(m — 1) 



(9) 



In the present work, we focus on the bipartite entangle- 
ment between the jth qubit and the remain k — 1 qubits 
constituting the fc-qubit system. Under the two mea- 



surements, it is denoted by Nj and Cj respectively. This 
corresponds to the case of m = 2 in (9), where the rela- 
tion reduces to Cj (p) > Afj (p) . And for the fc-qubit pure 
states Cj(\ip)) = J\fj(\ip)), which can be obtained directly 
from their definitions in (5) and (6). 

The entanglement invariants [39-41] of three-qubit 
pure states are five linear independent polynomials, 
which is invariance under LU transformations. For the 
three-qubit pure state p = \tp}(ip\, they are given by 



1, 
Z 5 



Trpf = l--rf, 



Ml 
2 

3Tr[( Pi ® Pj ) Pij ] - Trpf 
r 2 = (A^-^-C^) 2 , 



Trpf, (10) 



where pi — Tr jkp (i,j,k = 1,2,3) are the one-particle 
density matrices, pij = TikP are the two-particle den- 
sity matrices, and A/j and C- L j are the negativity and 
concurrence of the reduce states with the corresponding 
subscripts. The invariants Ii,2,3 are equivalent to the 
one-qubit linear entropies, which characterize the entan- 
glement between one qubit and the the remaining two 
qubits. The last invariant X 5 is equivalent to the three- 
tangle r, which describes the whole entanglement of the 
three-qubit system [42, 43]. The remaining invariant X4 
is related to the relative entropy of the two-qubit state 
Pij relative to the product state pi ® pj , and is a second 
measure of the entanglement of in the reduced states of 
Pij , independent of the concurrence Cij which can be de- 
termined by A/1,2,3 an d r [39]. 



C. SoDE and perturbation approach 

Under the local noise channels, the SoDE we focus on 
in the present article is the one between the ith qubit 
and the other part of the whole system. Utilizing the 
entanglement measure negativity, it can be expressed as 



Vi(p) 



dMi{p s ) 



dt 



(11) 



t=o 



where p s 



=1 £jP is the final state with p being the 



initial one. To compare with the widely utilized definition 
of the robustness by the ESD noise parameter, we present 
an alternative definition of robustness of entanglement, 
which is 



72, 4 = 1-Exp(-I7), 



(12) 



where T* — — is a characteristic time of the disentan- 
glement. 

We find that the SoDE 77, can be derived with the 
perturbation theory [44]. We called the method as the 
perturbation approach. In quantum mechanics, the per- 
turbation theory is applied to the systems whose Hamil- 
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tonian can be divided into 

H = Hq + eW, 



(13) 



where e is the smallness parameter, and H is the easily 
solvable unperturbed Hamiltonian. The eigenvalues and 
eigenstates of Hq are given by 



H \4> k ) = E^\4> k ). 



(14) 



The eigenvalues of the Hamiltonian H can be expanded 
in powers of the perturbation parameter e, 



E k = E^ + eE^' + e^> + ... 



(!) , ,2 P (2) 



(15) 



For the state |0fc) without degeneracy, the first approxi- 
mation is given by 



e£ ] - (4>k\w\<t> k ). 



(16) 



When the energy level is c^-fold degenerate, with 
the eigenstates |<^), (j = l,2...d k ), the corresponding 
eigenvalues of H in the first approximation are 



E k ,a = 4 0) + tPk,a, a -1,2.. .4, 



(17) 



where fik,a is the a-th eigenvalues of the dk x dk matrix 
W fe with the elements W k ,ij = {4> k \W\<f? k ) . One can no- 
tice that, when d k = 1, /3 k ,i becomes the result for the 
non-degenerate case in Eq. (16). And, the trace of W k 
satisfies 



TrW k = J^p k , a = J^(4\W\4). 

3 = 1 



(18) 



a=l 



Through the noise channel in an infinitesimal time dt, 
the state of a system can be written as 



p s = p + adt, 



(19) 



where p is the initial state, a is an Hcrmitian operator 
with Trer = 0. The negativity of the final state is given 
by jVi(ps) = J\Ti(p) — rjidt. Considering the partial trans- 
posed state pj* = p Ti + a Ti dt as the Hamiltonian of a 
quantum system, and dt as the smallness parameter, we 
find the SoDE r\i is determined by the first-order response 
of the eigenvalues of p Ti under the perturbation a Ti dt. 
Since the negativity is defined as the sum of the nega- 
tive eigenvalues of the partial transposed density matrix, 
only the eigenvectors of p Ti with the zero and the nega- 
tive eigenvalues contribute to the speed rji. Namely, the 
negativity of p s is 

M(Ps) = -2£ - 2 ]T X i- dt - 2 E A $,-^(20) 



where Xj t - are the negative eigenvalues of p Ti and X^__ 
are the first-order response under the perturbation a Ti dt, 



and A^q_ are the negative responses of the zero eigen- 
value. Therefore, the speed can be divided into two parts 
as 



(-) (o) 



(21) 



where r?| } = 2£.A^1 and r,f = -2£, Ag,. 
cording to the relations in Eqs. (16), (17) and (18) the 
first term can be expressed as 



,(°) 



Ac- 



,(-) - 



(22) 



where \ip k ), (k = l,2...d_) are the eigenvectors of p Ti 
with negative eigenvalues. And the second one corre- 
sponds to the eigenvectors \ipj) (j = 1,2. ..do) with the 



zero eigenvalues is 77-°' = J2i l^ol ~ > wnere A; q 1 

are the first-order responses of the zero eigenvalues. It is 
equivalent to 



(i) 



V, 



(0) 



Ti I 



Tra T % 



(23) 



where is the d x d matrix with the elements (^ mn = 



(demand || 
by ||G|| = TrVGGt. 



stands for the trace norm defined 



III. TWO-QUBIT ARBITRARY STATES 

First we explore the connection between the SoDE of 
a two-qubit quantum system and its entanglement prop- 
erties. Since the entanglement is invariant under permu- 
tations of the two qubits, we omit the subscripts of the 
entanglement and the SoDE. 



A. Pure states 

The pure states of a two-qubit system are always equiv- 
alent to 

\ip(6)) =cos0|OO) +sin0|ll), <9e[0,7r/ 4 ], (24) 

under LU transformations, with the entanglement 
jV[\ip{9))} = C[|V>(#))] = sin 29. The evolution of the neg- 
ativity for the states under the local depolarizing chan- 
nel can be calculated directly. But we are willing to give 
the the following steps to demonstrate the perturbation 
approach introduced in the above section. The partial 
transposed density matrix of the pure state (24) is given 
by 



cos 2 9 

















sin 26 
2 








sin 20 
2 

















sin 2 



(25) 
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which has no zero eigenvalue, and a eigenvector with the 
negative eigenvalue being \<f>~) — (|01) — 110))/\/2. Under 
the depolarizing channel in a finite time t, the partial 
transposed density matrix becomes 



l+s 2 +2scos28 
4 







4 

; sin 26 





'sin 20 

4 









l+s 2 -2scos 20 



,(26) 



where s = e *. The perturbation matrix can be derived 

as 



dt 



t=o 



cos 






2 9 





1/2 
sin 29 






sin 29 
1/2 








- sin 2 9 



(27) 



Substituting \<f> ) and a T into (21) and (22), we obtain 
the SoDE for the two-qubit pure states as 



t? = 2A/"+1, 



(28) 



which is a linear function of the entanglement. It is in- 
teresting to notice r\ — > 1 when Af — > 0, which indicates 
the ESD occurs in the pure state of this model even with 
a slight entanglement. In such a situation, the ESD time 
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equals to the characteristic time Tesd 
the robustness in (12) is given by IZ^ - 
the result given in [19]. 
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FIG. 1: Plots of the relations between the SoDE and neg- 
ativity of randomly generated two-qubit states: (a) 30000 
randomly generated states and the curves showing the up- 
per bound (solid) and the lower bound (dashed) of r/; (b) 
5000 randomly generated states weighted with the curve of 
the lower bound. 



For the arbitrary two-qubit states, since their entan- 
glement can't be completely described by one entangle- 
ment measure, it is very difficult to derive an analytical 
expression for the relation between the SoDE and the en- 
tanglement properties. Therefore we explore the bounds 
of SoDE for a given value of negativity and the states 
achieve the bounds, which is called as the frontier states 
in this article. 

Following the approach in [45, 46], we generate ran- 
domly a great number of two-qubit states and plot them 
in the rj-Af plane as shown in Fig. 1 (a). We fortunately 
find the pure states are the frontier states with the maxi- 
mum of SoDE. To derive the frontier states achieving the 
lower bound, we begin with the ansatz states [45, 46] 



Pansatz 



X + 






i 

a 

6 

y+i 



(29) 



where x, y, a, b, 7 > and x + y + a + b + -f = 1. The 
region of the ansatz states in the rj-Af plane is the same 
as the one of the arbitrary states. By trying to adjust the 
constrains imposed on the ansatz states without leaving 



the lower bound, we find the frontier states p m on the 
lower bound when x = y = b = 0, namely 

Pm = 7 |V>(7r/4))<V(V4)| + (1 - 7)|01)(01|, (30) 

with the pure states |^(7r/4)) defined in (24). To verify 
the frontier state p m , we randomly generate weighted 
random states, namely the mixtures of random states 
and p m with random weights, and plot them in the ij-Af 
plane. As shown in Fig. 1 (b), the region of physically 
acceptable states is encircled perfectly by the curve of 
the state p m . The results of the two families of frontier 
states presented the lower and upper bounds of the SoDE 
in two-qubit system as 



Af 2 + ^2Af{Af+l) 
1 + 2Af - ^2Af{N + 1) 



< 77 < 2Af+ 1. 



(31) 



When Af — > 0, the minimum of the SoDE approaches 
V2Af, and the ESD time and the characteristic time 
Tesd -> T* -> ^W/2. 

According with the results in [47], the two families of 
frontier states exactly is the ones in the comparison of the 
concurrence and negativity. Specifically, the pure states 
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FIG. 2: Plots of 30000 randomly generated states in (a) the 
£i — Xi plane and (b) the £2 — X2 plane. The solid curves are 
the results of p m in (a) and pk in (b). 



which has the maximal SoDE achieves the minimal con- 
currence when the negativity is fixed, and the states p m 
with the minimum of SoDE has the maximum of concur- 
rence for a given negativity. This results indicate that the 
concurrence reduces the attenuation of the negativity. 

Then, a question arises: Is there a definite link be- 
tween the extreme cases in the two topics? To explore 
the relations, we introduce four non-negative quantities 
as a = C - M, xi = 2A/" + 1 - V, 62 = M - M mm and 
X2 = V — rj mm , where -q mm is the lower bound of the 
SoDE in (31) and M mm = y/C 2 + (1 -Cf+C- 1 is the 
minimal negativity for a fixed concurrence satisfied by 
the frontier states p m . Plotting the randomly generated 
states in the £i~Xi and £,2~X2 planes as shown in Fig. 2, in 
the same approach for r\- Af, we obtain the frontier states 
in the two planes. They are the states p m for the first 
case, and 

Pk = -y|^(7r/4)><V(7r/4)| + o(|01><01| + |10)<10|), (32) 

with 7 + 2a = 1, for the second, which are shown by 
the solid lines in Fig. 2. By straightforward calculations, 
one can find that, £i(p m ) is a single- valued function of 
Xi(Pm), and equally true for £ 2 (Pfc) and X2(Pfe)- Hence, 
the states p m has the maximums of £1 for a given xi- 
When Xi(Pm) = 0, it can be derived that £i(p m ) = 0. 



Therefore, for arbitrary states, xi = =>■ £1 = 0. In 
other words, the concurrence reaching the minimum is a 
necessary condition for SoDE with the maximum. Simi- 
larly, one can find that the concurrence reaching the max- 
imum is a necessary condition for SoDE with the mini- 
mum. These demonstrates that there are other factors 
that affect the SoDE in two-qubit system in addition to 
the concurrence and negativity. 



C. Two-parameter states 

To study the influence on the SoDE by the concurrence 
and other quantities related with entanglement, we con- 
sider three classes of two-parameter mixed states in the 
similar form as the ansatz states (29). 

As an analytical example to show the relation between 
the SoDE and the concurrence, the first two-parameter 
class of states we studied is 

PC = 7|^/4))<V(7r/4)| + o|01)<01| + 6|10><10j, (33) 

with 7,8,6 e [0, 1] and 7 + a + b = 1, which is the ansatz 
state with x = y = 0. When a = b, it returns to the 
states pk in (32), whose concurrence and negativity have 
the same values. When ab — 0, it reduces the states with 
the maximal concurrence, such as p m in (30). Its SoDE 
can be expressed in terms of Af and C as 

„-2Af+i 2(1-C)(C-AQ(1+AQ 

which can returns to the bounds in (31) in the two cases 
of a = b and ab = respectively. In the Fig. 3 (a), 
one can notice in this family of states, when the other 
entanglement measure is fixed, the SoDE increases as the 
negativity increases, but decrease with the concurrence 
increasing. 

Considering the frontier states p m in (30) are precisely 
the maximally entangled mixed states in the region of 
C > 2/3 [45], we attempt to find the relation between 
the SoDE and the the degree of mixture, measured by 
the linear entropy as Sl(p) = f (1 — Trp 2 ). Since, for a 
given concurrence the pure states has the minimal linear 
entropy, we choose the hybrid states of (24) and (30) as 



Ps 1 =7lW>W)l+(l-7)|01)<01| 



(35) 



In Fig. 3 (b), one can notice for a fixed mixture in these 
states, the relation between the SoDE and negativity is 
approximately linearity. However, the SoDE does not 
clearly relate with the linear entropy S^. 

Rather than the mixture of the whole state, we are 
interested in the quantities describing the relations be- 
tween the subsystems. The total correlation in a bipar- 
tite quantum system has been defined as the difference 
between the sum of the von Neumann entropies of the 
two subsystems and that of the whole system, called the 
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FIG. 3: Plots of the relations between the SoDE and the 
negativity of the states (a) pc, (b) ps L and (c) pi tot , with the 
contours of the concurrence C, the linear entropy Sl and the 
quantum mutual information Ttot respectively. 



two-parameter states as 

Phot = Pansatz\y=Q,C=f(N')i 



(37) 



where f(x) = —x/2 + wx + 5x 2 /4. Here, actually, the 
choice of the function f(x) is quite arbitrary only requires 
f(x) > x when x £ [0, 1]. With the above form of the 
function, the region of the states pj tot in (37) in the -q-N 
plane is large enough, and the influence of the quantum 
mutual information can be shown clearly. In Fig. 3 (c), 
one can notice for a given Af, of course the concurrence 
C is fixed, the SoDE decreases with the Itot increasing. 
That is, the correlations other than the entanglement 
described by the concurrence also reduce the SoDE under 
a noise environment. 



IV. THREE-QUBIT PURE STATES 

In this section, we explore the SoDE of the three-qubit 
pure states, whose entanglement properties can be de- 
scribed by the invariants in (10). For a three-qubit sys- 
tem, there are two inequivalent classes of genuine tripar- 
tite entanglement, which are the GHZ state and the W 
state 

1 



|GffZ> = - ? =(|000) + |lll», 



|wO = -^=(|ooi) 



|010) + |100». 



(38) 



They can't be transformed each other by stochastic local 
operations and classical communication (SLOCC) [48]. 
The GHZ state possesses a maximum tripartite entan- 
glement characterized by the three-tangle [42] , for which 
case t{\GHZ)) = 1, but the residual bipartite entangle- 
ment of GHZ state is zero. However, the W state pos- 
sesses zero three-tangle, with r(|VF)) = 0, and maximizes 
the residual bipartite entanglement. 

In [19], the authors find the most robust symmetrical 
three-qubit pure states under the partially depolarizing 
channel are the GHZ-like states [see below Eq. (40)], and 
ascribe the robustness to the three-tangle. Our question 
is whether the conclusion is established when we con- 
sider the SoDE as the signature for the stability of the 
entanglement. Furthermore, what role do the five entan- 
glement invariants play in the SoDE respectively? 



quantum mutual information [2] 

lt ot {p)=S(p A )+S(p B )-S{p), (36) 

where S(p) — — Tr(crlncr) and /?^(and pb) is the reduced 
density operator. To explore the influence by the corre- 
lations outside of the entanglement, we choose a class of 



A. Symmetrical states 

To answer the first question, we start from the symmet- 
rical pure states, which are invariable under the permu- 
tations of the three particles. In this part, the subscripts 
of the negativity and the SoDE are also omitted. Gener- 
ally, an arbitrary symmetrical three-qubit entangled pure 
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(a) 




(b) 



FIG. 4: Plot of (a) 30000 randomly generated three-qubit 
symmetrical pure states; (b) 30000 randomly generated three- 
qubit arbitrary pure states in the r\ — M plane. The curves 
show the upper bounds (solid) and the lower bounds (dashed) 
in the two cases. 



the main result of the present subsection given in Eq. 
(48). We remark that the states | J) are not equal to the 
W-like states in [19] under LU transformation. 



For the state of \G), the partial transpose over the 
first qubit is p T G = a|000)(000| + — a)(|011><100| + 

1 100) <011|) + (1 - o)|lll)(lll|. The eigenstate with the 
negative eigenvalue can be obtained as \4>q) = (| 100} — 
|011))/\/2, and the ones with zero eigenvalue are |001), 
|010), 1 101) and 1 1 10) . The corresponding matrix in 
(23) is positive semidefinite, which is written in the basis 
{|001),|010),|101),|110)} as 



a G,0 ~ 



- 

2 u 

- 

u 2 











— 

2 U 





(41) 



Therefore, the eigenvectors with the zero eigenvalue do 
not have contribution to the SoDE, r](°\\G)) = 0. Simi- 
larly, ?/ )(| J)) = 0. Their SoDE are only determined by 
the responses of the negative eigenvalues under the noise 
channels, and can be derived as 



v(\G)) = ur+-, 



ri(\J)) = -Af+l, 



(42) 



which are the lower and upper bounds respectively. 
Hence, in the viewpoint of the SoDE, the GHZ-type 
states are more robust than the other symmetrical states, 
which fits with the result given in [19]. 



state can be written as [19] 

|$) =t 1 \000)+t 2 \W)+t 3 \W')+U\lll), (39) 

where ti (i = 1, 2, 3, 4) is the complex constant, J2i l^| 2 = 
1, and \W) = (a x <g> <7 X <g> cr^lVF) with a x being the first 
Pauli operator. The first three entanglement invariants 
satisfy X\ = I2 = X3 = 1 — , where N is the negativity 
between one qubit and its complementary two-qubit sub- 
system. Only three invariants in the symmetrical states 
are independent. 

As the two-qubit arbitrary states case, we can deter- 
mined the bounds of the SoDE in the present case. The 
region of the symmetrical three-qubit pure states in the 
rj- Af plane is shown by the randomly generated states in 
Fig. 4 (a). Following the approach in the case of the arbi- 
trary two-qubit states, we obtain the two frontier states 
on the lower and upper bounds are the GHZ-type states 
and the W-type states respectively 



\G) = v^lOOO) + vT^o|lll), 

I j) = Vb\w) + VT^ 



(40) 



However, it is very attractive to us that whether the 
robustness comes from the three-tangle. We can obtain 
the invariants in the two families of states as 



ATQG)) =2y/(l -a)a, ATQJ)) = -V2y/b(3 - 2b), 
t(\G))= 4(1 -a)a, r(|J)) = — = , (43) 



3^3 ' 

Ii{\G)) = 1 - 3a + 3a 2 , I 4 (| J)) = 1 - 3b + 4b 2 - 



16b 3 



For a given value of negativity, they satisfy t(\G)) > 
t(|J)) and Ii(\G)) > 1^(1 J)), in which the equalities 
hold when N = or N = 1. This suggests that there 
exists another candidate, the fourth invariant Z4 , for the 
main factor affecting the stability of the negativity. 



To answer the above question, we consider the hybrid 
states of \G) and |J), as 



|T) = Cl |000)+c 2 |W)+c 3 |lll), 



(44) 



where a, b e [0, 1]. They can be proved strictly by using 



where Cj (j — 1,2,3) are real and ^ - c 2 = 1. Its inde- 
pendent entanglement invariants (regardless of the nor- 
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malization condition) can be derived as 



General states 



TV 



6c 2 . c 2 , 



i 4 = 4 + H4 + ^ + 44 + 4 + \44{^2 + V3c 3 ), 



c 3 (4V3cl + 9clc 3 ) 



(45) 



Transposing the states |T)(T| partially on the first qubit, 
we derive the unique eigenvector with the negative eigen- 
value as 



37V + 4cl 



+ 



6cic 3 
37V + 6c? 



1 000) 



C2 



2v/3. 



cic 2 



A/3C3 

2 |011) — 



|001) 



C2 



V3c. 



C3 



|010) 



44 V3Af 



2V3c 



C2C3 



■| 100) 



3AT+ 24 + 6c| |101) _ 3AA + 2c|+6c| |no) 



6c 3 



6cic 3 



+ |111), 



(46) 



which is non-normalized. Under the depolarizing channel 
in an infinitesimal time, the form of a T can be derived 
directly. Substituting it and the eigenvector into (21), 
we obtain the SoDE in terms of the coefficients Cj as 
v(°ij c 2, C3). In this case, ?/ ) = means the eigenvectors 
with the zero eigenvalue do not have contribution to the 
SoDE. Considering the relations r = TV 2 and I4 = 1 — 
3/V 2 /4, satisfied by the GHZ- type states, we suppose the 
SoDE of the states |T) can be written as 

t ] = 3M+^ + (t-M 2 )X+(1 4 -1 + jM 2 )Y, (47) 

where X and Y are two undetermined functions of TV, t 
and I4. Substituting the invariants (45) into the above 
expression, and comparing it with 77(01,02,03), we obtain 
the form of SoDE for the state |T) as 



32 - 32I 4 - 127V 2 + 847V 3 + 697V 4 + 3r 2 
247V 2 (7V+ 1) 



■ (48) 



Very fortunately, we find the relation is also satisfied by 
the general symmetry states (39). 



From the above relation between the SoDE and the 
entanglement invariants, one can clearly notices that, r\ 
decreases with the increase of I4, but increases as the 
three-tangle increasing. In other words, the entangle- 
ment among the three qubits, quantified by the three- 
tangle, can enhances the SoDE, but the two-qubit corre- 
lations described by X 4 reduces the speed of negativity. 
Therefore, it is I4 but not the three-tangle r making the 
GHZ-type states to be the most robust symmetric pure 
states. In the analysis of the Ref. [19], the states they 
considered have no more than two independent parame- 
ters, which muddles up the influences of Z 4 and r. 



The discovery of the analytic expression of the rela- 
tion between the SoDE and the entanglement invariants 
in symmetrical states in (48) reveals the SoDE is an ef- 
fective tool to explore the role of different entanglement 
in the stability of multipartite entanglement. This moti- 
vates us to extend it to more general case, and study the 
influences of the other invariants on the SoDE. 

The general form of three-qubit pure states is given by 

I*) = co|000) + Cl |001) + c 2 |010) + c 3 |011) 

+c 4 |100) + c 5 |101) + c 6 |110) + c 7 |lll), (49) 

where Ci(i — 0,1, 2,..., 7) are complex constants and 
J2i \ c i\ 2 = 1- Without loss of generality, the analysis can 
be restricted to the negativity between the first qubit and 
the subsystem containing the other two qubits. 

Following our schedule in the symmetrical states case, 
we first derive the bounds of the SoDE for the arbitrary 
pure states. Plotting the randomly generated three-qubit 
pure states in the 771 -TV plane, we find only the lower 
bound are different with the one of the symmetrical states 
in the region of TVi > 1/2. And the altered lower bound 
corresponds to nothing but the result of the pure two- 
qubit state in (28). Thus, the frontier states on the lower 
bound of the SoDE for three-qubit pure states are the 
GHZ-type states when TV < 1/2, and the direct product 
states of an entangled state containing the first qubit 
and a single partite state of the remaining qubit, such as 
|/i) = |^(0))|O) with the |V>(6>)) defined in (24). Besides 
the symmetric pure states | J), we also find the frontier 
states on the upper bound contain the family of states as 



|n) = a|O)|V(0o))+/3|l)|V'(0i)), 



(50) 



where \a\ 2 + \/3\ 2 = 1, = (I ® a x )\i>{9)) and (0„ - 

tt/4)(0i - tt/4) = 0. 

To realize the ultimate aim of the SoDE for the arbi- 
trary pure states, we consider the subclass of states of ^ 
as 



|A) =c |000) + ci|001) + c 6 |110), 



(51) 



where Cj are restricted in the real numbers. It can be 
regarded a hybrid state of the two frontier states on the 
lower bound of the SoDE, which becomes the GHZ-type 
state when c = and returns to \fi) with c\ = 0. Its five 
entanglement invariants can be derived as 



x 4 

T 



x 2 = 1 - ^ 



1 - 24 4 



TV 2 
2 



6> 
„4 2 



(52) 



Cq 3cqC-^ -\- 3cqC-^ + c-y -\- Cg, 



4c? 4, 



where TVl = 2\J (c§ + 4)4- ^ n * ne ^ ne perturbation ap- 
proach, the first term of the SoDE for |A) the can be 
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gotten as 

-S-'(|A»=|V, + 1- "-^ 1 -^ . (53) 

There exist four eigenstates of the partial trans- 
posed density matrix of |A) with the zero eigenvalue, 
which span the degenerate subspace { 1 1 1 1) , c x 1 100) — 
co|101),|011),ci|000)-co|001)}. And the matrix cr^ 1 de- 
fined in (23) in the subspace is given by 














CqC 6 





(c|+cf) 


V c o+ c ? 


C C6 

V c o + c ? 




ol 
2 






2 



(54) 



The criterion for the matrix a 1 having a negative eigen- 
value is cfc 2 , — Cqc| < 0. From the relations in (52), we 



find it can be replaced by = (X 2 — X 3 ) 2 — r 2 /4 > 0, 
which is verified as the universal criterion for nonzero 
Tii for arbitrary three-qubit pure states. By using 
the relations in (10), the criterion can be rewritten as 
|ri2 — T13I — r > 0, where = Cfj are the two-tangle. 
The second part of the SoDE for |A) can be written as 

[ 0, 6 < , 

m (l A )) = I V(i-i3) 2 +J/?e-(i-i 3 ) Q>0 (55) 

On the other hand, we calculate a class of states with- 
out the three-tangle as 

= c |000) + ci|001) + C2IOIO) + c 4 |100), (56) 

where^ ■ | Cj | 2 — 1 with j = 0,1,2,4. It can be con- 
sidered as a generalization of the W-like states in [19] 
Similarly, for following the perturbation approach 
and replacing cj by the invariants, we get the SoDE 



v -16 - 32X 4 + 36A/" 2 + 847V-? + 57Af? - 12(X 2 - X 3 ) 2 + 12(X 2 + X 3 )(2 - Aff 
m ^ = 24A/?(M + 1) ; 



y/M2 + (J 2 -J 3 W-M 

where M = ±(5 - 3X X - 3X 2 - 3X 3 + 4X 4 ) = g( 4 - 6X 2 - 6X3 + 8X4 + 37V] 2 ). 

Based on an overall consideration of the analytic results of SoDE in the states |<&), |A) and we guess the form 
of SoDE for general three-qubit pure states to be 

( v -16 - 129 - 32X4 + 36A/ -2 + 84^ + 57A/? + 12(X 2 + X 3 )(2 - 7V" 2 ) 
m " } > ~ 247V 2 (7V 1 + 1) 

r e<o ; 

- < y/Mt+ATfe-M (58) 

I A?? ' « > • 

I 



By the numerical validation, the above unified form al- 
most is fulfilled by most all the three-qubit pure states. 
Namely, we take Aij = 1 771 — r) m \, where r\ m is the nu- 
merical solution by the finite difference method with 
At = 1Q- 9 . We generate 200000 sets of data for the At), 
which show that A77 is less than 10~ 5 when TVi > 10 -5 . 
When Afi < 10~ 5 , since the emergence of the phe- 
nomenon of ESD, the finite difference method is no longer 
suitable. 

From the form in (58), one can notice the conclusions 
about the influences of X 4 and r on the SoDE drawn from 
the symmetrical states also exist in the general case. The 
sum of the invariants X 2 and X 3 can enhance the SoDE. In 
the other words, the decay of the negativity between the 
first qubit and the subsystem containing the other two 
qubits, can be restrained by the other two negativity with 
the coequal sense. The unbalance between the second 



and the third qubits, shown by the differences (X 2 — X 3 ) 2 
and \ti2 — ti 3 |, can reduce the SoDE of the first qubit. 
Especially when |n 2 — ti 3 | > r, the vanishing negativity 
can be partially brought back, which reflected in nonzero 

(o) 

m ■ 



V. MULTIQUBIT STATES 

For the multiqubit system, in which the SoDE will be 
more complex, due to the absence of a uniform measure 
of entanglement, we only study some special symmetrical 
states. 

The fc-qubit GHZ-type states have been investigated 
widely as a family of exemplary states in the topic about 
the decay of entanglement in multiparty systems. It is 
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given by 



|G) fe =a|0)® fe +/3|l) 



(59) 



where a, (3 are complex constant and |a| 2 + |/3| 2 = 1. For 
the case of k > 3, the SoDE under the local depolarizing 
channel can be calculated directly as 



r,(\G) k ) = kAf- 



1 

2' 



(60) 



where the negativity J\f — 2\af3\. And, the alternative 
robustness is given by 



TZ v {\G) k ) = l-Exp( - 



2N 



2kAf+l 



(61) 



It is interesting to notice the SoDE given by (28) of the 
two-qubit GHZ-type states isn't contained in the form of 
(60). Based on the results in (42), we present a possible 
explanation as, the slop in the SoDE (28) of the two- 
qubit state comes from its similarity with the GHZ-typc 
states, but the intercept is related with its property of the 
W-type states, since they are the two different generaliza- 
tions of the two-qubit pure states to the three-qubit sys- 
tem. In the limit of large number of particles k — > +oo, 
the robustness (61) and the corresponding characteristic 
time approach TZ V — > 1/k and T* — > 1/k, which con- 
forms to the time at which such entanglement becomes 
arbitrarily small [16]. 

The second class of fc-qubit entangled states we con- 
sidered are the W states 



1 



|W) fc = -^(|00...01) 



|00...10) 



|10...00)), (62) 



whose negativity is J\f = 2\fk 
SoDE can be obtained as 



1/k. The corresponding 



= (k + 2)Vk—l + 2(k l)-(k-2)Vk^2 
k 

It is plotted in Fig. 5 (a), one can find that the SoDE 
r\ increases with k for 2 < k < 4, but decreases with k 
when k > 5. However, the robustness of \W)k only de- 
crease with k. When the number of particles k — > +oo, 
the negativity M — > 2/Vk, and the SoDE approaches a 
constant i](\W)k) —> 2. And then the alternative robust- 
ness of \W) k is given by K v (\W) k ) = 1 - Exp(-J\f/rj) -> 
1/Vk. It is worth mentioning the result in [49] that, for 
the fc-qubit system initial with the W states under the 
amplitude-damping channel, the negativity of the least 
balanced partitions decays with 1 / yk. 

When the number of qubits k — 3, the GHZ states 
\GHZ) k = (\0}® k + |l)® fc )/V2 and the W states in (62) 
become the two frontier states (40) for the symmetric 
case with a = 1/2 and 6 = 1. And the negativity of their 
superposition 



(a) 




\Z(q,<p))k = y/q\GHZ) k -e i *y/T=q\W) k , (64) 



FIG. 5: (a) Plot of the relation between the SoDE and the 
number of qubit for \W) k - (b) Plot of the maximal influence 
of ip on the SoDE of the states \Z(q,<fi)) k for k — 3 (solid), 
k — 4 (dashed) and k = 5 (dot-dashed). 



does not depend on the phase factor ip [19]. Therefore, 
the influence of ip on the SoDE is related with the en- 
tanglement properties described by r and Z4. Since, for 
arbitrary values of k, the negativity of \Z(q, <p>)) k 

., V(" fc2 + 6fc " 4 )« 2 + 2 ( fc2 - 5fc + 4 )9 + 4(fc - 1) 
Al = .(65) 

also only depends on the parameter q, we can consider 
the maximal effect of <p 

Ar,(q) = max{r)[\Z{q, cp a )} k ] - r)[\Z(q, >p b )) k ]}, (66) 

as a characteristic quantity to indicate the affection of 
the entanglement besides the negativity. In Fig. 5 (b), 
we plot the relation between A.r](q) and q for k = 3, 4, 5. 
It is shown that, when k — 3, the phase tp presents a 
significant affection on the SoDE for arbitrary q. But, for 
the case of k = 4, the influence is not distinct when q < 
0.4, where the proportion of the GHZ state in \Z(q, (f))k 
is less than the one of the W state. When k = 5, the 
values of Ai](q) can hardly be seen in the Fig. 5 (b). With 
a numerical simulation, we find the amount of Ar)(q) is 
no more than 10 -4 , for k = 5 and 6. That is, the SoDE of 
the states \Z(q, ip))k are almost entirely on the negativity 
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when k becomes large. 

By contrast, one can choose the GHZ-type states (59) 
with the negativity N = 2^/k — 1/fc, and derive the dif- 
ference between the values of SoDE in (60) and (63). It is 
easy to find the difference increases with increasing num- 
ber of qubits, when k > 4. This reveals, the affection 
of the entanglement besides the negativity on the SoDE 
becomes more prominent as the number k increases. Ac- 
cording to these results, the family of the states \Z(q, ip))k 
is not a good sample to explore the role of different entan- 
glement components in the stability of the entanglement 
in the systems with a large number of particles. 



VI. CONCLUSION AND DISCUSSION 

In the present study, we investigate the dynamical 
properties of entanglement of multiqubit systems under 
local partially depolarizing channels. In this model, each 
qubit is independent couples with its own environment. 
Our main concern is the bipartite entanglement between 
one qubit and its complementary subsystem, measured 
by negativity. The relations between the SoDE and the 
entanglement properties in an entangled state are ex- 
plored. 

In two-qubit system, for the pure states, we get the 
analytical expression of the SoDE, which is determined 
completely by the negativity. For the arbitrary states, 
using the form of the ansatz states, we gain the upper 
and lower bounds of the SoDE. The pure states has the 
maximal SoDE, and the mixture of mutually orthogonal a 
Bell state and a separable pure state as (30) achieves the 
minimum of SoDE. The minimum and the maximum of 
concurrence are shown to be the necessary conditions for 
the two bounds respectively. With the aid of some classes 
of two-parameter states, we find the SoDE can be reduced 
by both the concurrence and the total correlation. 

In three-qubit system, we derive the analytical expres- 
sion of the SoDE in terms of the entanglement invariants 
for arbitrary pure states. The GHZ-type states are shown 
to be the most robust, which is consistent with the result 
in [19]. However, by the relation between the SoDE and 
the invariants, we find the main reason for the robustness 
in the GHZ-type states is the two-qubit correlation de- 
scribed by the fourth invariant I4. On the contrary, the 
three-tangle among the whole system can enhance the 
speed of negativity under the local depolarizing chan- 
nel. The decay of the negativity can be restrained by 
the other two negativity with the coequal sense. The 
unbalance between two qubits can reduce the SoDE of 
the remaining one in the three-qubit system. Especially 
when |ti 2 — T13I > r, the vanishing negativity of the first 
qubit can be partially brought back. 

At last, we study fc-qubit system with the aid of some 
families of exemplary states. For the fc-qubit systems in 
an arbitrary superposition of GHZ state and W state, 
the influence of the entanglement outside the negativity 
becomes less noticeable when k increases to five or six. 



In the limit of large number of particles, for the GHZ- 
type states, the characteristic time and the correspond- 
ing robustness defined in (12) base on the SoDE, is in- 
versely proportional to k. This coincide with the result of 
the time at which such entanglement becomes arbitrarily 
small given in [16]. Under the same condition, the robust- 
ness of the fe-qubit W states approaches 1/y/k. A similar 
behavior has been reported in [49], exhibited by the the 
negativity of the least balanced partitions, in the fc-qubit 
system initial with the W states under the amplitude- 
damping channel. These results show the effectiveness of 
the SoDE and the corresponding alternative robustness 
to quantize the stability of the entanglement in multipar- 
tite quantum systems. 

Finally, we briefly discuss the universality of our def- 
initions and results, especially about the choices of the 
entanglement measures and the noise channels. The main 
reason for adopting the negativity as the entanglement 
measure in our present work about the SoDE is its com- 
putability, and fatherly its speed being also computable 
effectively with the aid of perturbation theory. A com- 
putable quantity always has more advantages for physi- 
cists. Because of the necessity and sufficiency of N > 
for entanglement in pure states, the SoDE defined in (11) 
and the related alternative robustness can be used to ex- 
plore the stability of the entanglement in arbitrary di- 
mensional multiparty pure states. In addition, although 
M > only is the sufficiency condition for entanglement 
in mixed states with the dimension larger than six, the 
speed of negativity can still be considered as a signa- 
ture for stability of the entanglement in the states with 
a nonzero Af. Meanwhile, the robustness related with 
the ESD noise parameter is suspect, because its practi- 
cal computation is usually based on the PPT criterion. 

In principle, on can choose any measures of entangle- 
ment and study the corresponding SoDE. We conjecture 
that, in the multiparty system with pure states, the qual- 
itative conclusions for other entanglement measure are 
similar with the ones for negativity. 

Taking the concurrence in (7) for instance, one can 
define its speed as the form in (11) 



From the relation in (9), it is directly to find the speed 
of the concurrence for a pure state in multiqubit system 
satisfies 

V C <V N , (68) 

where rj N denotes the corresponding SoDE in (11). Es- 
pecially, for the two-qubit case, 

if = 2C + 1 = V N . (69) 
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And, for the fc-qubit GHZ-type states 

fcC + | <v C <kC + ± = V N , (70) 

where the upper bound comes from (68), and the lower 
bound can be derived as following. After passage through 
the depolarizing channel, the final state can be divided 
as 

pf = [1-S(s)}p e (s)+S(s)p d (s), (71) 

where p e (s) is a state in the space of {|0),|1)} <g> 
{lO)®^- 1 ), (l)®^- 1 )}, p d (s) is a separable diagonal state, 
and they together with the non-negative coefficient S(s) 
depend on the noise parameter s. By the convexity of 
concurrence, 

C(pf)<[l-S(s)]C[p e (s)} (72) 

which leads to the lower bound in (70). Moreover, based 
on the bounds in (70), we can obtain the corresponding 
alternative robustness ~R£ = 1 — Exp(— Cjrf) — > 1/fc, 
when fc — > oo, which is the same as the result for nega- 
tivity. Then, the above results left a interesting question 
for us: Can we find a explicitly expression for rf of the 
multiqubit pure states? 

However, for the mixed states, the speed of concur- 
rence has a problem of singularity. For example, in the 
two-qubit system, the form of ansatz states (29) is held 
under the local depolarizing channel. The speed of con- 



currence can be obtained as rp = y £ || + vTlt _ ft > 
which approaches infinity when a > and b — > (or 
a — > and b > 0). This property also exists under other 
local operations [50] , which holds the form of the ansatz 
states and drastically enforces full rank of the states when 
a = or b = 0. However the corresponding alternative 
robustness TZ£ has a finite value. 

In order to compare with the depolarizing channel, we 
also calculate the SoDE of the GHZ-type states \G)h 
under the local dephasing channel, which is desc ribed 
by the Kraus operators as Eq = y(l + e~ t )/2I and 
Ei = - e~ t )/2a 3 . The SoDE for both the concur- 
rence and the negativity are 

v = kC = kAf, (73) 

which also lead the alternative robustness lZ n — > 1/fc, 
with fc — »■ oo. This is in accord with the conclusion about 
the limit of large number of particles in [25] . 



Acknowledgments 

We are grateful to the reviewer for valuable comments. 
F.L.Z. thanks Jing-Ling Chen, Fu-Guo Deng and Bao- 
Kui Zhao for their valuable discussions and encourage- 
ment. This work is supported by NSF of China (Grant 
No. 11105097). 



[1] A. Einstein, B. Podolsky, and N. Rosen, Phys. Rev. 47, 
777 (1935). 

[2] M. Nielsen and I. Chuang, Quantum computation and in- 
formation (Cambridge University Press, Cambridge, UK, 
2000). 

[3] A. Ekert, Phys. Rev. Lett. 67, 661 (1991). 

[4] C. Bennett, G. Brassard, and N. Mermin, Phys. Rev. 

Lett. 68, 557 (1992). 
[5] C. Bennett, G. Brassard, C. Crepeau, R. Jozsa, A. Peres, 

and W. Wootters, Phys. Rev. Lett. 70, 1895 (1993). 
[6] C. H. Bennett and S. J. Wiesner, Phys. Rev. Lett. 69, 

2881 (1992). 

[7] T. Yu and J. H. Eberly, Phys. Rev. Lett. 93, 140404 
(2004). 

[8] M. Yonac, T. Yu, and J. Eberly, J. Phys. B: At., Mol. 

Opt. Phys. 39, S621 (2006). 
[9] A. Al-Qasimi and D. James, Phys. Rev. A 77, 012117 

(2008). 

[10] T. Yu and J. Eberly, Science 323, 598 (2009). 

[11] K. Roszak, P. Horodecki, and R. Horodecki, Phys. Rev. 

A 81, 042308 (2010). 
[12] J. Laurat, K. Choi, H. Deng, C. Chou, and H. Kimble, 

Phys. Rev. Lett. 99, 180504 (2007). 
[13] M. Almeida, F. De Melo, M. Hor-Meyll, A. Salles, S. Wal- 

born, P. Ribeiro, and L. Davidovich, Science 316, 579 

(2007). 

[14] A. Salles, F. de Melo, M. P. Almeida, M. Hor-Meyll, S. P. 



Walborn, P. H. Souto Ribeiro, and L. Davidovich, Phys. 

Rev. A 78, 022322 (2008). 
[15] C. Simon and J. Kempe, Phys. Rev. A 65, 052327 (2002). 
[16] L. Aolita, R. Chaves, D. Cavalcanti, A. Acin, and 

L. Davidovich, Phys. Rev. Lett. 100, 80501 (2008). 
[17] Z. Man, Y. Xia, and N. An, Phys. Rev. A 78, 064301 

(2008). 

[18] Z. Liu and H. Fan, Phys. Rev. A 79, 064305 (2009). 
[19] B. Zhao and F. Deng, Phys. Rev. A 82, 014301 (2010). 
[20] K. Zyczkowski, P. Horodecki, M. Horodecki, and 

R. Horodecki, Phys. Rev. A 65, 012101 (2001). 
[21] G. Vidal and R. Tarrach, Phys. Rev. A 59, 141 (1999). 
[22] Y. Sagi, I. Almog, and N. Davidson, Phys. Rev. Lett. 

105, 053201 (2010). 
[23] J. Du, X. Rong, N. Zhao, Y. Wang, J. Yang, and R. Liu, 

Nature 461, 1265 (2009). 
[24] T. Cubitt, F. Verstraete, and J. Cirac, Physical Review 

A 71, 052308 (2005). 
[25] M. Hein, W. Diir, and H.-J. Briegel, Phys. Rev. A 71, 

032350 (2005). 

[26] G. Vidal and R. F. Werner, Phys. Rev. A 65, 032314 
(2002). 

[27] A. Peres, Phys. Rev. Lett. 77, 1413 (1996). 

[28] K. Zyczkowski, P. Horodecki, A. Sanpera, and M. Lewen- 

stein, Phys. Rev. A 58, 883 (1998). 
[29] A. Borras, A. P. Majtey, A. R. Plastino, M. Casas, and 

A. Plastino, Phys. Rev. A 79, 022108 (2009). 



14 



[30] Y. S. Weinstein, Phys. Rev. A 82, 032326 (2010). 

[31] W. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 

[32] K. Chen, S. Albeverio, and S.-M. Fei, Phys. Rev. Lett. 
95, 040504 (2005). 

[33] Kraus, States, effects, and operations: fundamental no- 
tions of quantum theory (Springer- Verlag, Berlin, 1983). 

[34] W. Diir and H. Briegel, Phys. Rev. Lett. 92, 180403 
(2004). 

[35] C. H. Bennett, D. P. DiVincenzo, J. A. Smolin, and W. K. 

Wootters, Phys. Rev. A 54, 3824 (1996). 
[36] C. H. Bennett, G. Brassard, S. Popescu, B. Schumacher, 

J. A. Smolin, and W. K. Wootters, Phys. Rev. Lett. 76, 

722 (1996). 

[37] V. Vedral, M. B. Plenio, K. Jacobs, and P. L. Knight, 

Phys. Rev. A 56, 4452 (1997). 
[38] M. Horodecki, P. Horodecki, and R. Horodecki, Phys. 

Lett. A 223, 1 (1996). 
[39] A. Sudbery, J. Phys. A: Math. Gen. 34, 643 (2001). 
[40] A. Acin, A. Andrianov, L. Costa, E. Jane, J. Latorre, 

and R. Tarrach, Phys. Rev. Lett. 85, 1560 (2000). 
[41] A. Acin, A. Andrianov, E. Jane, and R. Tarrach, J. Phys. 



A: Math. Gen. 34, 6725 (2001). 
[42] V. Coffman, J. Kundu, and W. Wootters, Phys. Rev. A 

61, 052306 (2000). 
[43] R. Lohmayer, A. Osterloh, J. Siewert, and A. Uhlmann, 

Phys. Rev. Lett. 97, 260502 (2006). 
[44] W. Greiner, Quantum Mechanics : An Introduction 

(Springer, 2000), p. 273, 4th ed. 
[45] W. J. Munro, D. F. V. James, A. G. White, and P. G. 

Kwiat, Phys. Rev. A 64, 030302 (2001). 
[46] T. Wei, K. Nemoto, P. Goldbart, P. Kwiat, W. Munro, 

and F. Verstraete, Phys. Rev. A 67, 022110 (2003). 
[47] F. Verstraete, K. Audenaert, J. Dehaene, and B. Moor, 

J. Phys. A: Math. Gen. 34, 10327 (2001). 
[48] W. Diir, G. Vidal, and J. Cirac, Phys. Rev. A 62, 062314 

(2000). 

[49] R. Chaves and L. Davidovich, Physical Review A 82, 
052308 (2010). 

[50] F. Verstraete, J. Dehaene, and B. DeMoor, Physical Re- 
view A 64, 010101 (2001). 



